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We develop a Hamiltonian formulation of Bianchi type-I cosmological model in 
conformal gravity, i.e., the theory described by the Lagrangian £ = C ' a bcdC abcd , 
which involves the quadratic curvature invariant constructed from the Weyl ten- 
sor, in a four-dimensional spacetime. We derive the explicit forms of the super- 
Hamiltonian and the constraint expressing the conformal invariance of the theory, 
and we write down the system of canonical equations. To seek out exact solu- 
tions to this system we add extra constraints on the canonical variables and we 
go through a global involution algorithm that possibly leads to the closure of the 
constraint algebra. This enables us to extract all possible particular solutions that 
may be written in closed analytical form. On the other hand, probing the local 
analytical structure we show that the system does not possess the Painleve prop- 
erty (presence of movable logarithms) and that it is therefore not integrable. We 
stress that there is a very fruitful interplay of local integrability- related methods 
such as the Painleve test and global techniques such as the involution algorithm. 
Strictly speaking, we demonstrate that the global involution algorithm has proven 
to be exhaustive in the search for exact solutions. The conformal relationship of 
the solutions, or absence thereof, with Einstein spaces is highlighted. 



1 Introduction 

Unlike what happens in general relativity, where unwanted second-order deriva- 
tive terms of the metric can be discarded from the gravitational action through 
pure divergences, there is no chance whatsoever to weed out such terms in the 
context of higher-order gravity theories. It is fortunate however that a consis- 
tent method of building up a Hamiltonian formulation of those theories does 
actually exist; it is a generalization of the classical Ostrogradsky formalism. El 
Basically it consists in introducing auxiliary degrees of freedom that encom- 
pass each of the successive derivative terms higher than first order. Without 
resorting explicitly to this method Boulware worked out a Hamiltonian formu- 
lation of quadratic gravity. □ Contradistinctively, Buchbinder and Lyakhovich 
developed a canonical formalism for the most general quadratic gravitational 
Lagrangian in four dimensions, by employing the aforementioned generalized 
Ostrogradsky method. I In a previous work we have already applied Boul- 
ware's canonical formalism to Bianchi cosmologies, for the pure R 2 variant 
of the general quadratic theory. □ In this contribution we consider the con- 
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formally invariant case, which is based on the quadratic Lagrangian density 
£ = C abcd C abcd , where C a bcd is the Weyl tensor.^ Making use of a slightly 
different generalized Ostrograsky construction (as compared to Buchbinder and 
Lyakhovich's formalism) we derive the canonical form of the conformally in- 
variant action, which was first obtained by Boulware. We further particularize 
the Hamiltonian formalism to the Bianchi type-I cosmological model and write 
down the canonical equations. This is our starting point for seeking out exact 
particular solutions and for asking whether the conformal Bianchi type-I cos- 
mological model is integrable or not. Jn that respect, we sum up very recent 
results which can be found elsewhere. 



2 Conformal gravity and Bach equations 

Consider a four-dimensional spacetime (A4,g) and the quadratic gravitational 
action 

s = -\ I d 4 xV=gC abcd c abcd , (l) 

4 JM 

where C a bcd is the Weyl tensor. Its variation with respect to the metric yields 
the conformally invariant fourth-order equations 

B a b '■= 2V m V™C ma (,„ + C ma bnR mn — 0, (2) 

which have been put up by BachB who adopted Weyl's paradigm of a con- 
formally invariant gravitational theory, El though without considering the addi- 
tional Weyl 1-form. In Eq. (||), B a b is called the Bach tensor; it is symmetric, 
trace-free and conformally invariant of weig ht -1. Although Eq. (§) is a very 
compact formula it is not appropriate for calculational purpose with computer 
algebra — the contracted double covariant derivative of the Weyl tensor is 
rather heavy to compute even for simple metrics. Recently Tsantilis et al have 
provided an algorithm for the MathTensor package that gives the Bach equa- 
tions in a much more tractable form, especially with regard to cosmological 
applications, and which is based on the decomposition of the Riemann tensor 
in its irreducible pieces. □ In accordance with their results we write down the 
following equivalent expression of B a b-, as given by Eq. (^|), 

B a b = — n (^Rab — ~^9ab^j + - V a Vfei£+ (C ma b n + R ma bn + Rmb9an) R™™ ■ (3) 

The simplest cosmological model exhibiting non trivial physical degrees of 
freedom in the conformally invariant gravitational theory based on the action 

a We adopt rWald's conventions of sign and definitions of curvature tensors (but the extrinsic 
cur vature) . cl 
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given in Eq. is the spatially homogeneous anisotropic Bianchi type-I model. 
(The isotropic flrw cosmological models are conformally flat.) Writing the 
metric in such a way that the conformal invariance becomes manifest already 
from the outset, 



(4) 



the corresponding Bach equations, as given in Eq. (||), can be derived with 
the help flf symbolic computational packages such as the excalc package in 

REDUCE, u 



3 Hamiltonian Bianchi type I cosmology 

3. 1 Hamiltonian formalism and canonical equations 

The conformally invariant action in Eq. (|TJ) can be cast into Hamiltonian form 
by means of a generalized Ostrogradsky construction. 

Assume first that spacetime is foliated into a family of Cauchy hyper- 
surfaces Tit of unit normal n a . The induced metric h a b onto these hyper- 
surfaces is defined by the formula h a b = g a b + n a rib. The way the hyper- 
surfaces are embedded into spacetime is provided by the extrinsic curvature 
tensor K a b :— —\£nh a b, where £n denotes the Lie derivative along the nor- 
mal n a (we adopt Cft as a generalized notion of time differentiation). The 
standard ADM variables are introduced: the lapse function N and shift vector 
A\ The 3 + 1 -splitting of spacetime enables us to express the Lagrangian 
density £ = — zV~9 C a bcdC abcd only in terms of the quantities defined onto 
St. After some algebra we obtain the following equations, with respect to the 
ADM basisfl 



ijk 
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(5) 
(6) 



Owing to the fact that the Weyl tensor identically vanishes in three dimensions, 
the Lagrangian density in Eq. (m) reduces to 



-NVh (2C ninj C ninj 



n ... r ni i k \ 



(7) 



'Latin indices i,j,k,... refer to spacelike components in Et; the ADM basis is defined by 
e n := n = (dt — N l di)/N and Si = di\ a stroke '|' denotes the covariant derivative onto St, 
and the symbol n indicates a component along e n . 
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where C n y& and C„i n j are given by Eq. (||) and Eq. (||), respectively. 

Now consider that the induced metric and the extrinsic curvature Kij 
are independent variables — i.e., are introduced as auxiliary Ostrogradsky 
variables. In order to recover the definition i*Cy = —^Cfihij we must trade the 
original Lagrangian density for a constrained Lagrangian density 

£ = N- X Z + A 4 -? {Cnhij + 2K. l0 ) , (8) 

with Lagrange multipliers \ lJ as additional pzariables. Thus we must resort 
to Dirac's formalism for constrained systems. E3 The generalized Ostrogradsky 
construction enables us to cast the action in Eq. into canonical form, 

S= f d 4 xN [p« Cfthij + Q iJ C H K i:i - S) c (h, K, p, Q)] , (9) 

with the conjugate momenta p^ — A y and Q y = —2^/hC nm \ and where the 
canonical Hamiltonian density is given by 

Sj c = -2p l U< t] +VhC ni0k C^ k -^^-Q l \ ij -Q^^R t3 -KK l] Q^. (10) 

Dirac's constraint analysis yields, besides the usual super-Hamiltonian and 
super-momentum constraints, one first-class constraint that is the generator of 
conformal transformations. It reads explicitly 

x = 2p + K kl Q kl »0. (11) 

Moreover we can get rid of the spurious canonical variables Q and K; only 
the traceless part of the corresponding tensors remain as relevant canonical 
variables. This is consistent with the fact that conformal gravity exhibits six 
degrees of freedom. 

Performing a canonical transformation that we have defined in a previous 
workta in order to disentangle terms stemming respectively from the pure R 2 
and conformal variants of the general quadratic theory, we get 

S= j d 4 x [lL> + U+$+ + n_/L + Q+V+ + Q-V- - NH C - A cV J , (12) 

where the first-class constraints TL C rs (super-Hamiltonian) and (p c ss are 
given respectively by the following expressions 

H c = — L [n+-p+ + Il-V- + 2Q+ (Vl - V 2 _) - 4V+V-Q-] 
V6 

-\e-^(Ql + Q 2 _), (13) 
< Pc =Tl f ,-V + Q + -V-Q-. (14) 
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A suitable gauge- fixing condition that eliminates variables /i and and the 
choice N — yield the final form of the canonical action 



s = [ d 4 x [n +/ 3+ + n_/L + q+t + + Q-V- - n c ] , (15) 

Jm 1 J 
where the super-Hamiltonian is now given by the following expression 

H c = 4= Q- - U+V+ - - 2Q+ {V% - V 2 )] - ]■ (Q\ + Q 2 _) . 

VD * 

(16) 

Varying the action as given by Eq. ( |l5| ) with respect to the remaining canonical 
variables and their conjugate momenta, we obtain the canonical equations for 
the Bianchi type-I model in conformal gravity ^] 

P± = -^±> ( 17 ) 

n± = o, (is) 

Q+ = -^(U + +4V-Q--4V + Q+), (19) 

Q- = --4= (II_ + 4"P_Q+ + AV+Q-) , (20) 
V6 

f + = 2-{Vl-Vl) + Q+, (21) 

V- = -^V+V- + Q- (22) 
V6 

Instead of the fourth-order Bach equations, we have now at our disposal the 
nice differential system, given by Eq. @ to Eq. (||), which is more appropri- 
ate for applying singularity analysis methods (Painleve test) in order to extract 
all the exact solutions. In contrast with these methods which probe the local 
analytical structure of the canonical system, we can also seek out exact solu- 
tions by performing a global involution algorithm on specific extra constraints 
chosen in accordance with local results from the analytic structure. 



3.2 Global involution of extra constraints 

The involution method consists in applying the Dirac-Bergmann consistency 
algorithm on our system, with the Poisson brackets defined with respect to the 

c To enable a straightforward comparison with Boulware's formalism we also have performed 
a canonical transformation that interchanges the coordinates and momenta V± and Q±. 
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canonical variables j3±, H±, Q±, V±, and after suitable conditions have been 
imposed. Strictly speaking, the steps of the global involution algorithm are 
the following: 

• Impose an appropriate extra constraint on the canonical variables; 

• Require that constraint to be preserved when time evolution is consid- 
ered. This gives rise to secondary constraints and possibly to the deter- 
mination of the Lagrange multiplier associated with the extra constraint; 

• Repeat the second step (involution) until no new information comes out. 

Once the involution algorithm has been performed we can classify all the con- 
straints into first class and second class and proceed further to the analysis of 
the particular system. 

Among the set of constraints we have considered in our analysis two are 
more significant. The first expresses that the ratio of the variables V± is 
constant. Any solution to the canonical equations that satisfies that specific 
constraint is conformally equivalent to an Einstein space. Consistency of the 
extra constraint yields only one secondary constraint and both are second 
class. We eliminate the associated spurious degrees of freedom and reduce the 

) to one binomial equation of Briot and Bouquet. 
The representations of the solution on thejjeal axis are complicated expressions 
involving the Weierstrass elliptic function.lj As a particular case of that analysis 
we obtain the general axisymmetric solution. The second interesting constraint 
requires that the momenta IT± be zero. In accordance with local results from 
the analytic structure, the general solution to the complete system with n± = 
is easy to produce under analytic form, upon integrating linear ODE's of the 
Lame type. It provides a specific example of a solution to the Bach equations 
that is not conformally related to an Einstein space. □ 

4 Analytic structure of Bianchi type I cosmology 

The existence of the above solutions, whether particular solutions of the gen- 
eral differential system or general solutions of specialized systems, tells nothing 
about the integrability or non-integrability of the complete system and gives 
no information whatsoever about the mere accessibility of an exact and closed- 
form analytic expression of its general solution. This is due to the fact that 
the global involution algorithm of the extra constraints, as operated above, is 
not related with integrability and may even prove to be non-exhaustive. We 
have tackled the integrability issue through an invariant investigation method 



canonical equations (jl£ 22 
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of intrinsic properties of the general solution.^] In particular, we have proved 
analytically that the system under consideration is not integrable: its general 
solution exhibits, in complex time, an infinite number of logarithmic transcen- 
dental essential movable singularities — i.e., an analytic structure not compat- 
ible with integrability in the practical sense; the quest for generic, exact and 
closed-form analytic expressions of the solution is hopeless. We stress that this 
result holds under spacetime transformations within the equivalence class of 
the Painleve property. Eil 



5 Conformal relationship with Einstein spaces 

The vanishing of the Bach tensor is a necessary condition for a Riemann space 
to be conformal to an Einstein space. It means that any solution of vacuum 
general relativity or any space conformal to an Einstein space are also a solu- 
tions of conformal gravity. The converse however is not true: there exist spaces 
not conformally related to Einstein spaces that satisfy the Bach equations. 113 
A necessary and sufficient condition for a space to be conformal to an Einstein 
space is the existence of a function <r(x), i.e., conformal factor, that satisfies 
the differential equations 

1 R 

Lab = V a (jV 6< T - V a V b< 7 - -g ab g cd Vc<jV d <J - ^ e2a 9ab = 0, (23) 

where the tensor L a b is defined by L a b := (Rg a b — 6i? a fc)/12. For the Bianchi 



type-I model it is not difficult to show that the above conditions in Eq. (23) 
imply precisely that the ratio of variables V± is constant (our first extra con- 
straint in the algorithm). In that case the conformal factor can be uniquely 
determined as a function involving the Weierstrass elliptic function. Moreover 
it is possible to obtain the explicit form of the constant scalar curvature, R, 
of the conformal Einstein space. On the other hand, our solution obtained 
by imposing H± = is not conformal to an Einstein space, for the ratio of 
variables V± is not constant in that case. We have thus confirmed explicitly 
Schmidt's conjecture_of the existence of what he calls 'non-trivial' solutions to 
the Bach equations. E3 
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